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Abstract. We study the blow-up criterion of smooth solutions to the 3D MHD equations. 
By means of the Littlewood-Paley decomposition, we prove a Beale-Kato-Majda type blow-up 
criterion of smooth solutions via the vorticity of velocity only, i. e. supj &z J Q T \\Aj(V x^Woodt, 
where A- is a frequency localization on |£| ~ 2 J . 

Key words. MHD equations, Blow-up criterion, Littlewood-Paley decomposition 

AMS subject classifications. 76W05 35B65 

1 Introduction 

We consider the 3D incompressible magneto-hydrodynamics (MHD) equations 

du . __ __ l__, l9 , __,, 

vAu + u-Vu = -Vp W + b ■ Vfe, 

dt F 2 ' 

db 

— — rjAb + u ■ Vb = b ■ Vu, (11) 

V -u = V-6 = 0, 

, u(0, x) = u (x), 6(0, x) = b (x). 



(MHD) 



Here u, b describe the flow velocity vector and the magnetic field vector respectively, p is 
a scalar pressure, v > is the kinematic viscosity and r] > is the magnetic diffusivity, 
while uo and bo are the given initial velocity and initial magnetic field respectively, with 
V • no = V • bo = 0. If v = 7] = 0, is called the ideal MHD equations. 

Using the standard energy method, it can be easily proved that for given initial data 
(uo,bo) € H S (M. 3 ) with s > ^, there exists a positive time T = T(\\(uo,bo)\\jj s ) and a unique 
smooth solution (u(t, x),b(t, x)) on [0,T) to the MHD equations satisfying 

(u, b) e C([0, T); H s ) n ^((0, T); H s ) n C((0, T)-H s+2 ). 

Whether smooth solutions of (jl.ip on [0, T) will lead to a singularity at t = T is an outstanding 
open problem, see Sermange and Temam |17j . Caflisch, Klapper and Steele [2] extended the 
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well-known result of Beale-Kato-Majda pQ for the incompressible Euler equations to the 3D 
ideal MHD equations, precisely, under the condition: 



Jo 



|V x u(t)||oo + ||V x 6(t)||oo)tft < oo, (1.2) 



then smooth solutions (u, b) can be extended beyond t = T. Recently, there are some re- 
searches which have refined (|1.2H such as 

/ (||Vx«(t)||^ + ||V X b(t)\\g )dt<oo, (see [21]). 

limsup / (||A i (Vx«)(t)[| 00 + [|A i (Vx6)(t)|| 00 )cft = <y<M J (see PI) 

jGZ JT-s 

for some positive constant M, and Aj is a frequency localization on |£| ~ 2 J . These results 
can be easily extended to (jl.ip with u, r\ > 0. Wu [22] also extended some Serrin type criterion 
for the Navier-Stokes equations to the MHD equations. Many relevant results can be found 
in [201 [2T] and references therein. 

However, some numerical experiments [TJ [16] seem to indicate that the velocity field plays 
the more important role than the magnetic field in the regularity theory of solutions to the 
MHD equations. Recently, inspired by Constantin and Fefferman initial work [6] where the 
regularity condition of the direction of vorticity was used to describe the regularity criterion 
to the Navier-Stokes equations, He and Xin [8] extended it to the MHD equations, but did 
not impose any condition on the magnetic field b which was consistent with the result of 
numerical experiments. Precisely, they showed that the solution remains smooth on [0,T] if 
the vorticity of the velocity w = V x u satisfies the following condition 

\w(x + y, t) — w(x, t)\ < K\w(x + y, t)\ \y\ 2 if \y\ < p \uu(x + y,t)\ > O, (1-3) 

for t € [0, T] and three positive constants K, p, f2. Also, they [8] and Zhou [25] obtained some 
integrability condition of the magnitude of the only velocity u alone, or the only gradient of 
the velocity Vu alone to characterize the regularity criterion to the MHD equations, i.e. 

f T 2 3 

/ \\u(t)\\ldt < 00, - + ~<1 3<p<oo; (1.4) 
Jo Q P 



or 

i-T 



f 2 3 3 

/ \\Vu(t)\\ldt < 00, - + -<2 -<p<oo. (1.5) 
Jo Q P 2 



We restrict ourselves to substitute Vu by the vorticity w in (jl.5p . In the case p < 00, using 
the Biot-Savart law [12] and the bounds of the Riesz transforms [18] on L p (l < p < 00), the 
condition (jl.5p can be replaced by 

r T 2 3 3 

/ IIV x u(t)\\ldt < 00, - + -<2, -<p<oo. (1.6) 
Jo Q P 2 
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However, since the lack of continuity of Riesz transforms on L°°, their results missed the 
important marginal case p = oo which exactly corresponds to the Beale-Kato-Majda criterion. 
In the case of the Euler equations, Beale, Kato, and Majda get around this difficulty by using 
the logarithmic Sobolev inequality: 

||Vw||oo < C(l + l|V x u|| 0O log(e+ |M|jf.))> s>5/2. (1.7) 

For a refined form of this inequality, it can be referred to I14j . In order to make use of (1.7), 
one need to estimate the higher order derivatives of the solution (at least in H s , s > 5/2). But, 
in the case of the MHD equations, it seems difficult to control them by the only ||V x u||oo. 
Therefore, as in [2], if the logarithmic Sobolev inequality (11. 7p is used, one can only derive a 
criterion described by the vorticity of u and b . This difficulty is avoided by the following two 
observations. On one hand, the H 1 norm of the solution can be used to control any H s norm 
of the solution, while the H 1 norm of the solution (u, b) can be controlled by ||V«||oo- On the 
other hand, if we make use of Littlewood-Paley decomposition to decompose the nonlinear 
terms into three parts: low frequency, middle frequency and high frequency, and deal with 
each part by virtue of different estimates, we can refine ||Vu||oo to ||V x ti||oo- It should 
be pointed out that we do not apply the Littlewood-Paley decomposition to the equation 
itself as some researches do before, since if we localize the equation on a dyadic partition, we 
cannot control the H 1 norm of the solution (u,b) via ||V x u\\loo in the end when summing 
up every dyadic partition. 

Finally, we remark that the blow-up criterion we will establish in the framework of mixed 
time-space Besov spaces may be the most relaxed in some sense as for the incompressible 
Euler equations [15] and the Idea MHD equations [3], where the losing estimate for the 
solution and the logarithmic Sobolev inequality are applied to set up the blow-up criterion, 
but in this paper, if we follow their method, as we mentioned above, we cannot characterize 
the blow-up of smooth solutions by V x it only. 

Now we state our result as follows. 

Theorem 1.1. Let (uo,&o) S H s , s > | with V • uq = V • 6o = 0. Suppose that (u,b) £ 
C([0,T);H s )r]C 1 ((0,T);H s )r]C((0,T);H s+2 ) is the smooth solution to OTTjh If there exists 
an absolute constant M > such that if 

limsup/ ||Aj(V x u)\\oodt = 5 < M, (1.8) 
jez Jt-e 

then 5 = 0, and the solution (u, b) can be extended past time t = T . In other words, if 

r-T 

fl 



limsup [ ||Aj(V x ujWoodt > M, 



then the solution blows up at t = T. Here Aj is a frequency localization on |£| ~ 2 3 , see 
Section 2. 

Remark 1.1. For the Navier-Stokes equations (with b = in U.l\) ). Kozono, Taniuchi fl 
and Kozono, Ogawa, Taniuchi 111]/ refined the Beale-Kato-Majda criterion to 

r-T r-T 

\V x u(t)\\BModt < oo and j ||V X it(t)||_g dt < oo 



o Jo 
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respectively. Here B^ ^ stands for the homogenous Besov spaces, see Section 2 for the 
definition. The condition U.8\) is weaker than all the above mentioned conditions. Hence, 
this also improves the results of JiOj, 1 1 If . For the further explanation of \ 1.8\) . it can be 
referred to the remarks after Theorem 1 in 115)/ . 



Remark 1.2. Very recently, Wu [23] use the energy estimate combined with the Bony para- 
product technique to derive many interesting regularity criterions for the generalized MHD 
equations. In the marginal case, the regularity criterion obtained there can be expressed as 

f T 

\\u(t)\\]^ S dt < oo, or / |lit(£)|l B i+e dt < oo, 

^OC,OC Jq ^00,00 

for some 5,e > 0. Here B pq stands for the inhomogenous Besov spaces, see Section 2 for the 
definition. But the velocity u cannot be replaced by its vorticity, since the Riesz transformation 
is not bounded in B^ x . 

Remark 1.3. For the Ideal MHD equations, whether similar result holds is still open, since 
the viscous term plays an important role in our proof. 

Notation: Throughout the paper, C stands for a "harmless" constant, and changes from 
line to line; || • \\ p denotes the norm of the Lebesgue space IP. 

2 Preliminaries 



Let <S(R ) be the Schwartz class of rapidly decreasing functions. Given / 6 <S(M ), its Fourier 
transform T f = / is defined by 
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/(0 = (2tt)-3 / e- lx <f(x)dx. 



Now let us recall the Littlewood-Paley decomposition (see [H [TU] ) . Choose two nonnegative 
radial functions x, (p S 5(M 3 ), supported respectively in B = {£ E M 3 , |£| < |} and C = {£ G 
R 3 , § < |f | < §} such that 

5>(2-'"0 = i, eei 3 \{0}. 

Let h = T~ x tp and h = T~ x x, the frequency localization operator is defined by 



A 3 f = (p(2-W)f = 2^ [ h{py)f{x - y)dy, 
S d f = X (2-W)f = 2^ [ ~h{Vy)f{x - y)dy. 



Informally, Aj is a frequency projection to the annulus {|£| ~ 2 3 }, while Sj is a frequency 
projection to the ball {|£| < 2-?}. Observe that A,- = Sj — Sj— i- Also, if / is an I? function 
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then Sjf — > in L 2 as j — > — oo and Sjf —> f in L 2 as j — ► +00 (this is an easy consequence 
of ParsevaPs theorem). By telescoping the series, we thus have the homogeneous Littlewood- 
Paley decomposition 

+00 

j=-oo 

for all / G L 2 , where the summation is in the L 2 sense. 

Let s G R, 1 < p, q < 00, the homogenous Besov space B^ q is defined by 



£p, 5 = {/ € ll/IU < 00} 



p.q 

where 

1 

a 



E 2iS9 ii A ;/n; 



-j=-oo 

(usual modification if q = 00), and Z'(M. 3 ) can be identified by the quotient space of S' /V 
with the polynomials space V . The inhomogenous Besov space Bp is defined by 

B s p>q = {/ G S'(R 3 ); ||/|| B . j? = ||5 /|| p + ||{2^|A,/|| P },>o||^ < 00}. 

We now denote the operator (I — A) 2 by A which is defined by 

A?(0 = (1 + lfl a )»/(fl. 
More generally, A s / for s£l can be identified with the Fourier Transform 

AV(0 = (i + l£lW(0- 

For s £R, we define 

and the Sobolev space iP is denoted by H s = {/ G 5'(IR 3 ); ||/||.H"« < °°}- The usual Sobolev 
space H s ' p is endowed with the norm 

||/|| H .., 4 ||a s /||l,. 

We can refer to |19| for more details. 



Lemma 2.1. Let k G N. There exist constants C independent of f , j such that for all 
1 < p < 9 < 00 

sup ||9 Q A,/|| g < C2 J ' fc+3j '^-i)||A i /|| p , (2.2) 

|a|=fc 

ll^/H, < C^H^Aj/Hp- (2-3) 
Here (k = 1, 2, 3) is i/ie Riesz transform in M 3 . 
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The proof of this lemma can be found in [U [13] . 

Remark 2.1. Suppose that the vector function f is divergence- free, and set g = V x f . Then 
there exist constants C independent of f such that 

||V/||p < C%|| p , VKp<oo. (2.4) 

If the frequency of f is restricted to some annulus {|£| ~ 2 J }, then there holds 

||V/|| P < C||<7|| p , Vl<p<oo. (2.5) 

Indeed, the inequality \2.J$ can be derived from the Biot-Savart law llStf and the bounds of 
the Riesz transforms Jigj/ on L p (l < p < oo), while the inequality 12. 5\) can be deduced from 
the Biot-Savart law and \2. 3\) . 

Lemma 2.2 (Commutator estimate). Let 1 < p < oo, s > 0. Assume that f,g£ H s ' p , then 
there exist constants C independent of f , g such that 

\\A s (fg)-fA s g\\ LP < C(|| V/||i^ |M| HS -i, P2 + \\f\\H^ 3 \\g\\ LPi ) (2.6) 

withp2,P3 G (l,+oo) such that 

111 11 

V Pi P2 P3 Pi 

This lemma is well-known and for a proof, see [9]. 

3 Proof of Theorem 1.1 

We will divide the proof of Theorem 1.1 into two steps. 

Step 1. H 1 estimates. 

In this step we will show there exists e > 0, 

sup {\\u(t)\\ m + \\b(t)\\ m ) < C{\\u(T - e)\\ m + \\b(T - e)\\ m + e). (3.1) 

te[T-e,T) 

Let w(t, x) = V x u(t, x) and J(t, x) = V x b(t, x). Taking the curl on both sides of (jl.ip . it 
can be written as 

dw 

— uAw + (n • V)w - (w • V)u - (b • V) J + (J • V)6 = 0, 

?r (3-2) 
— - nAJ + (u ■ V) J - (J • V)u - (b ■ V)w + (w ■ V)6 = 2T(b, u) 

with 

(d 2 b ■ d 3 u - d 3 b ■ d 2 u\ 
d 3 b ■ d lU - dib d 3 u\. 
d\b ■ d 2 u — d 2 b ■ d\u I 
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Multiplying the first equation of (|3.2[) by w, the second one of (13. 2|) by J, then adding the 
resulting equations yields that 

^(IH*)ll! + 11^)111) +^l|V^(i)||| + r ? ||VJ(i)|||= / (wV)u-wdx+ [ (J-V)u-Jdx 
2 at Js.3 J R 3 

-[ ((J -V)b-w + (wV)b- J)dx + 2 [ T(b,u) ■ Jdx 
= I + II + III + 2 IV, (3.3) 
where we have used the facts 

/ (u ■ V)w ■ wdx = / (u • V) J • Jdx = 

and 

((b ■ V) J • w + (6 • V)w ■ J)dx = 0, 



which can be deduced from the divu = div6 = and integrating by part. 

In what follows, we will deal with each term on the right hand side of (|3.3|) separately 
below. Let us begin with estimating the term /. Using the Littlewood-Paley decomposition 
(|2.ip . we have 

w = V] Ajto = V] Aj-u) + V] Ajw + ^ Aj-u;, (3.4) 

j€Z J<-JV -N<j<N j>N 

where N is a positive integer to be determined later. Putting (|3.4|) into / produces that 

N 

I = y / (w ■ V)it • Ajwdx + y / (w ■ V)u ■ Ajwdx + > (w ■ V)u ■ Ajwdx 



j<-N K j=-N K j>A f * 

I + I2 + /3. 



Using the Holder inequality, (|2.4p and (|2.2p to obtain that 

< IMMlVulla J] l|Ai^||oo < Clklll J] ^IIAHb < C^^lkHl, ( 3 - 5 ) 

j<—N j<-N 

\h\ < \\w\\2\\Vu\\2 ^2 HAjiolloo < C||w||| ^2 \\^j w \\oo- (3.6) 

-N<j<N -N<j<N 

From the Holder inequality, (|2.4p . (|2.2p and Gagliardo-Nirenberg inequality, it follows that 
\h\ < \\w\\ 3 \\Vu\\ 3 \\ A j w h < C\\w\\l J2 2 ' W A M\2 

j>N j>N 
1 1 

<C\\w\\l(Y, 2 ~ j Y(Y; 22j \\ A i w ^) 2 <C2~t\\ w \\ 2 \\Vw\\1 (3.7) 
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By summing up (I3.5p - (|3.7|) . we get 



|/| < C(2~~^ N \\w\\l + \\w\\l \\ A j w \\oo + 2~^\\w\\ 2 \\Vw\\l). (3.8) 

-N<j<N 



Using the Littlewood-Paley decomposition (|2.ip to V«, II can be written as 

TV 

11= V / (J- V)Aju -Jdx+y" [ (J- V)A jU -Jclx+S" [ (J- V)A jU ■ Jdx. 



j<-N 



j=-N 



j>N- 



Then the Holder inequality, (|2,2p , (|2,5p and Gagliardo-Nirenberg inequality allow us to show 
that 

\II\ <C(2-^ N \\w\\ 2 \\J\\l + \\J\\l ^2 \\Ajw\loo + 2~t[|J|| 2 ||Vw|| 2 ||VJ||2). (3.9) 

-N<j<N 

Similarly, using the Littlewood-Paley decomposition (12. ip . Ill and IV can be written re- 
spectively as 



111= V f J ■ (V6 + (V6) T ) • Ajwdx + V / J • (V6 + (V6) T ) • A jW dx 

j<-N Jk3 j=-N J ^ 

+ V / J • (V6 + (V6) T ) • Ajwdx, 

j>N Jr3 

N 

IV = V / T(6, A^u) • Jdx + V / T(6, Aj-u) • Jdx + V / T(6, Aj-u) • Jdx. 
Then exactly as in the derivation of (|3.8p . (|3.9p . we can deduce that 



\III\+2\IV\<C(2-^ N \\w\\ 2 \\J\\l + \\J\\l W A M\oc 

-N<j<N 

+2-t||j|| 2 ||Vw|| 2 ||VJ||2). 
Combining (l3T8|) - (l3TT0|) with (J33|), Young inequality yields that for t £ [0,T) 



(3.10) 



d_ 
~dt 



(imt)||I + ||J(t)|||)+2H|V^)||| + 2r ? ||VJ(i)|| 



< c 2-f^(iKt)||I + ||j(t)||I) + \\ A Mt)\U\\w(t)\\ 2 2 + \\m\i) 

-N<j<N 

N 



+ 2-~ (\\w(t)\\ 2 + || J(t)|| 2 ) (||V^(t)||| + ||VJ(t)||i; 



(3.11) 



Now let us choose a fixed positive integer N such that C2~t (||tt;(t)||2 + || J(t)|| 2 ) < min(z/, rj), 
i.e. 



N > 



■log 



C 



log 2 \min(z/, rf) 



■(iKt)ii2 + iiJ(t)ii 2 ; 



+ i. 



(3.12) 



where log + x = log(e + x). Thus (ETTT1) and (IXT2I) imply that for t 6 [0, T) 
~(||u;(t)||jj + II J(t)lll) + ^l|Vu;(t)||l + i7||VJ(t)||l 
<C^ \\AMt)\U\\w(t)g + \\J(t)\\ 2 2 )+C, (3.13) 

j=~N 

which together with the Gronwall inequality gives that for t G [0, T) 

\\w(t)\\ 2 + \\J(t)\\ 2 < exp (c V T ||A i u;(t / )||oo^ / V^+ Il^(°)ll2 + ll J ( )ll2)- 

V J= _ N Jo J 

Recalling the choice of N in (|3,12p . it follows from the above estimate that 

\\w(t)\\ 2 + ||J(i)|| 2 < exp fciog+ {\\w{t)\\ 2 + \\J(t)\\ 2 ) sup f HA^^OHoodt' 



x (VCt + H0)|| 2 + || J(0)|| 2 ). 
For simplicity, let E(t) = ||u>(i)|| 2 + || J(i)||2, the above inequality implies that 

sup Eit) < exp fciog+ ( sup E{t)) sup / ||A i u;(t / )[| 0O dt /> ) (VCT + E(p)). (3.14) 
[o,T) V [o,T) jez Jo / 

We point out that the inequality (|3.14p still holds if the time interval is replaced by [T — e, T). 
It follows from (3.14) that 

sup E(t) < exp I log + sup E(t) sup / || Aj-u^t'^loocft' ) 

te[T-e,T) \ te[T-s,T) jeZ JT-e ) 

x (VC~£+\\w{T-e)\\ 2 + \\J{T-e)\\ 2 ). 
Defining Z(T) = log ( sup E(t) + e), the above estimate means that 

te[T-e,T) 

Z{T) < log (VCs + E(T - e) + e) + CZ(T) sup / ||A^u;(t')|| 0O £ft / . (3.15) 

j€Z JT-e 

If we choose M = in Theorem 1.1, the condition (|1.8f) ensures that there exists a small 
positive £o such that 

f T 1 

Csup / ||Aj«;(t')||oorft' — > Ve e (0,e ). 
jez Jt- £ 2 

Then the inequality (|3.15p implies that 

Z(T) <2log(E(T-e)+e), VeG(0,e )- (3.16) 



On the other hand, it is easy to prove that the solution (u, b) satisfies the energy inequality 



«(f)ll! + Wmili + 2 / (^livu(t')lll + vW)\\ 2 2)dt' < M*)\\l + \\b(s)\\l 



which together with (|3.16p . (|2.4[) implies (|3.ip . 
Step 2. H s (s > 1) estimates. 

For completeness, we will show how to deduce H s estimates from H 1 estimates. Taking 
the operation A s on both sides of (jl.ip . multiplying (A s u, A s b) to the resulting equation, and 
integrating over R 3 , we get 

^|(l|A'«(*)lll + l|A s 6(t)||I) + i/||VA'u(*)ll2 + r/||VA s 6(t)||i 

= -/ A s {u-Vu)A s udx+ ( A s (b-Vb)A s udx- [ A s (u ■ Vb)A s bdx 

JR3 J R 3 J R 3 

+ / A s (fr- Vu)A s 6cfe. (3.17) 



Noting that divu = divb = and integrating by parts, we rewrite (|3.17p as 

^(l|AM*)lll + \m(t)\\t) + *l|VA'u(t)lll + »7||VA'6(t)||l 

= - f (A s (u-Vu)-u- A s Vu)A s udx - [ (A s (u • V6) - « • A s Vb)A s bdx 

+ / (A s (b ■ V6) - 6 • A s V6)A s n + (A s (& • Vu) - 6 • A s Vu)A s 
Jm. 3 



A 



? bdx 

ni + n 2 + n 3 . (3.18) 



For the term IIi, it follows from (|2.6jl . Holder inequality, Gagliardo-Nirenberg inequality and 
Young inequality that 

I Hx I < C(||V«||2[|Vtt[|#«-M + ||ii||^s,4||Vn|| 2 )||u||^ s ,4 

< (7||V«|| 2 ||«||| s ||Vn||| s < C||Vm||^||«||^ s + ^-\\Vu\\ 2 H s. (3.19) 

The other terms can be treated in the same way: 

|n 2 | + |n 3 | < C(||Vti|| 2 ||V6|| ff .-M + ||«|| h m||V6|| 2 )||6||h.,4 

+ C(||V6|| 2 [|V6[| h .-m + IHIhm||V6|| 2 )|M|^,4 

+ C(||V6|| 2 ||Vu|| H .-i,4 + ||6|| h .,4||V«|| 2 )||6||h.,4 
< C(||V^||1 + \\Vb\\l)(\\uf Ha + \\bf Hs ) + V -\\Vuf Ha + l\\Vbf Hs . (3.20) 

By summing up (|3.19p and (|3.20p with (|3.18p . we get 

j t (\Ht)\\h + \\m\\h)+Hv<t)\\h +ri\mt)f HS 



<c(l|Vu(t)||! + ||v&(t)|H)(||u(t)||| s + "--' 2 
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Then the Gronwall inequality yields that 

\\<t)\\h + \\b(t)\\h + /Vl|Vu(f)|&. + V \\Vb(t')f HS )dt' 
Jo 

<C(||u(0)||k + ||6(0)||^.) expfi snp \\(u(t'),b(t'))\\%\ (3.21) 

V t'€[0,t) / 

Hence we have the H s regularity for the solution at t = T and the solution can be continued 
after t = T. This completes the proof of Theorem 1.1. 
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